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CONTROLLABILITY OF RANDOM SYSTEMS: UNIVERSALITY 
AND MINIMAL CONTROLLABILITY 

SEAN O’ROURKE AND BEHROUZ TOURI 


Abstract. For a large class of random matrices A and vectors 6, we show that 
linear systems formed from the pair (A, b) are controllable with high proba¬ 
bility. Despite the fact that minimal controllability problems are, in general, 
NP-hard, we establish universality results for the minimal controllability of 
random systems. Our proof relies on the recent developments of Nguyen-Tao- 
Vu |12| concerning gaps between eigenvalues of random matrices. 


1. Introduction 

In recent years, the topic of controllability of linear systems has regained a great 
deal of attention due to its potential impact on the controllability of networked sys¬ 
tems. Several researchers have provided conditions on controllability of specific lin¬ 
ear systems with a given set of parameters. These works include [D El El uni HH US] 
where several necessary and sufficient conditions have been obtained for controlla¬ 
bility of a system with a given transition matrix and a given input vector/matrix. 
Also, it was shown in |14j that the problem of minimal controllability of a lin¬ 
ear system (i.e. finding the sparsest input vector that makes a given linear system 
controllable) is NP-hard, in general. 

Motivated by these works, we study the controllability of random systems. That 
is, we consider linear systems whose parameters (as far as controllability is con¬ 
cerned) are random. We confirm a common belief that “most systems are highly 
controllable” even when one deals with systems of a very discrete nature. 

For the purposes of this note, we define controllability in terms of Kalman’s rank 
condition [9]. 

Definition 1 (Controllable). Let A be an n x n matrix, and let & be a vector in 
M". We say the pair (A, b) is controllable if the n x n matrix 

{b Ab ■■■ A"-i&) (1) 

has full rank (that is, rank n). Here, the matrix in (El is the matrix with columns 
b, A6,...,A"-i6. 

Our work is an attempt to address the following commonly believed conjecture. 

Conjecture 1. As n —> oo, the probability that (A, e^) is controllable for all i G 
{!,...,n} approaches one, where A is the adjacency matrix of an Erdos-Renyi 
random graph G{n,p), and Ci is the ith member of the standard basis in M". 

In addition, we also discuss the following conjecture from [Tj^ . 

^The conjecture in [7] is slightly different than the one stated here as it involves the uniform 
probability distribution on graphs with n vertices. 


1 



2 


SEAN O’ROURKE AND BEHROUZ TOURI 


Conjecture 2 ([!]). As n ^ oo, the probability that {A, !„) is controllable ap¬ 
proaches one, where A is the adjacency matrix of an Erdos-Renyi random graph 
G{n,p) and 1„ is the all-one vector in R.". 

We prove a stronger version of the first conjecture and provide some insights 
into the second conjecture. Our work heavily relies on the recent advancements 
in inverse Littlewood-Offord theory (see |13) for a survey) and its implications on 
eigenvalue spacing of random matrices recently studied by Nguyen, Tao, and Vu 
[12) . We prove that, under a very general assumption on the random matrix A 
and the (deterministic or random) vector b, the pair {A, b) is controllable with high 
probability. We also establish that, for a general class of random matrices A, the 
pair {A, Ci) is controllable for all i £ {1,... , n} with high probability. Our results 
suggest that, perhaps, the problem of controllability and minimal controllability is 
not a concerning problem for large systems. 

1.1. Outline and notation. The paper is organized as follows. We introduce the 
main objects of study in Section |21 and we state our main results in Section [31 We 
then present the technical tools we will use to prove our main results in Section [d] 
The proofs of our main results are presented in SectionjS) Finally, Section|6]contains 
a few closing remarks and some related open problems. 

Throughout the paper, we use the following notation. For any matrix A, ||j4|| 
is the spectral norm of A. A matrix obtained by deleting the ith row and the zth 
column of A is called a first minor of A. Since we only work with first minors, we 
simply refer to this matrix as the ith minor of A. We let A"^ be the transpose of 
A, and A* is the conjugate transpose of A. 

For a,n n X n Hermitian matrix A, we let 

< • ■ ■ ^ A„(A) 

denote the (ordered) eigenvalues of A (counting algebraic multiplicity) with corre¬ 
sponding orthonormal eigenvectors ui(A),. .. ,Vn{A). Since the eigenvectors of A 
are not uniquely determined, unless otherwise stated, we always assume that the 
collection {ui(A),... ,'?;„(A)} is any orthonormal basis of eigenvectors of A such 
that 

Avi{A) = \i{A)vi{A). 

We emphasis the fact that, throughout this paper, all the eigenvectors under con¬ 
sideration are unit length. 

We let ei,..., e„ denote the n elements of the standard basis in R". If u and v 
are vectors in C" (alternatively R"), then ||u|| is the Euclidean norm of u, and u*v 
(respectively u^v) denotes the standard Euclidean inner product of u and v. 

For an event E, we let E'^ denote the complement of E. For a random variable 
f with a given distribution, if there is no confusion, we refer to the distribution of 
^ simply as f. 

Throughout this paper, we consider controllability of single-input systems. We 
often refer to a linear system whose transition matrix A or input vector b are 
random, as a random system. 

2. Random matrix ensembles 

The focus of this paper is on controllability of a random matrix A and a (deter¬ 
ministic or random) vector b. In particular, we address the following question. 
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Question 1. Let b S K." be nonzero. If A is an n x n random Hermitian matrix, 
what is the probability that {A,b) is eontrollable? 

Of course, the answer to Question [T] depends on the distribution of A and the 
particular choice of b. In this note, we consider two important ensembles of real 
symmetric random matrices. 

2.1. Wigner random matrices. Wigner [19] introduced the following class of 
real symmetric random matrices with independent entries in the 1950s. 

Definition 2 (Wigner matrix). Let ^ and C be real random variables. We say 
W = is an n X n Wigner matrix with atom variables ^ and C if hb is a 

real symmetric matrix whose entries satisfy the following: 

• the entries {wij '■ 1 < i < j < n} are independent random variables, 

• the upper triangular entries {wij : 1 < i < j < n} are iid (independent and 
identically distributed) copies of 

• the diagonal entries {wu : 1 < i < n} are iid copies of f. 

A classical example of a Wigner matrix is the adjacency matrix of a random 
graph. Let G(n,p) be the Erdds-Renyi random graph on n vertices with edge 
density p, i.e. it is a graph on n vertices such that each pair of vertices {i,j} is 
an edge in G(ji,p) with probability p, independent of other edges. We denote by 
A = the zero-one adjacency matrix of G{n,p). That is. 



In particular, the upper-diagonal entries of A are iid copies of a zero-one Bernoulli 
random variable (with parameter p), and the diagonal entries are all zero. 

In general, we will be interested in Wigner random matrices whose entries are 
sub-gaussian. 

Definition 3 (Sub-gaussian). We say the random variable ^ is sub-gaussian with 
sub-gaussian moment k > 0 if 


IPdCI >t)<K ^exp(—Kt^) 


for all t > 0. 

There are many examples of sub-gaussian random variables. For instance, Gauss¬ 
ian random variables are sub-gaussian. In addition, all bounded random vari¬ 
ables (i.e. random variables that are supported on a bounded subset of R) are 
sub-gaussian. 

2.2. Orthogonally invariant ensembles. We also consider a class of orthogo¬ 
nally invariant ensembles. 

Definition 4 (Orthogonally invariant ensemble). We say the nxn random symmet¬ 
ric matrix M is drawn from an orthogonally invariant ensemble if M has probability 
distribution 


f‘{dM) = -;^exp(—trQ(M))c?M 


(2) 
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on the space of n x n real symmetric matrices, where dM refers to the Lebesgue 
measure on the n(n +1)/2 algebraically independent entries of M, Q is a real-valued 
function, and is the normalization constant 




e-trQWrfM. 


Here Q{M) is defined by the functional calculus for M. 


The prototypical example of an orthogonally invariant ensemble is the Gaussian 
orthogonal ensemble (GOE). The GOE is defined by the probability distribution 
in ([2]) when Q{x) = . In this case, a matrix drawn from the GOE is actually 

a Wigner random matrix. Indeed, if M = {rnij)2j=i is drawn from the GOE, the 
elements {iriij '■ 1 < i < j < n,} are independent Gaussian random variables with 
mean zero and variance 1 -I- 6ij , where 


<5 




1 , 

0 , if i j. 


We refer the reader to [3] for further examples and properties of orthogonally in¬ 
variant ensembles. 


3. Main results 

In this section, we present the main results of this work. To improve readability, 
we present the results first, and we postpone the technical details and proofs until 
the subsequent sections. Our main results below address Question [T] for the random 
matrix ensembles introduced above. In particular, we will consider two distinct 
cases: when b is deterministic and when 6 is a random vector. 

3.1. Results for Wigner matrices and deterministic vectors. We first con¬ 
sider the controllability of (IE, e^) when IE is a Wigner matrix. 

Theorem 1 (Minimal controllability of Wigner matrices). Let ^ be a sub-gaussian 
random variable with mean zero and unit variance; let C, be a (possibly degenerate) 
sub-gaussian random variable with mean zero. Assume W is an n x n lEi^ner 
random matrix with atom variables ^ and Let F be an n x n deterministie real 
symmetric matrix such that ||E|| < n'^ for some fixed parameter 7 > 0. Then, 
for any a > 0, there exists C > 0 (depending only on a, 7 , and () such that 
(W F, a) is controllable for every 1 < i <n with probability at least 1 — Cn~°‘. 

Remark 1. The matrix M ■.= W -\- F can be viewed as a Wigner matrix whose 
entries are not required to have mean zero. Alternatively, it is often useful to view 
M as a random perturbation of the matrix F. 

As a consequence of Theorem [U we obtain the following corollary for the adja¬ 
cency matrix of a random graph. 

Corollary 1 (Minimal controllability of random graphs). Fix 0 < p < 1 and a > 0. 
Let A be the adjaceney matrix of G{n,p). Then there exists C > 0 (depending only 
on p and a) such that {A, Ci) is controllable for each 1 < i < n with probability at 
least 1 — Cn~°‘. 

As a result of this corollary, Conjecture [T] follows immediately. This result also 
suggests that for large systems, the difficult problem of minimal controllability m 
is most probably a trivial problem. 
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3.2. Results for Wigner matrices and random vectors. We now consider 
the case when 6 is a random vector. In particular, we show below (see Corollary[2]) 
that {A, h) is controllable with high probability when A is the adjacency matrix of 
a random graph G{n,p) and & is a random binary vector. 

We begin, as before, with a result for a very general class of Wigner matrices. 
Indeed, Theorem [2] below shows that, with high probability, (W, &) is controllable 
when IT is a Wigner matrix and 6 is a random vector, independent of IT. More 
generally, we also consider deterministic perturbations of IT and b. That is, we 
show (IT + F,b + ii) is controllable with high probability, where F is a deterministic 
matrix and /r is a deterministic vector. 

Theorem 2. Let ^ be a sub-gaussian random variable with mean zero and unit 
variance; let C, be a (possibly degenerate) sub-gaussian random variable with mean 
zero. Assume W is an n x n Wigner random matrix with atom variables ^ and 

Let F be an n X n deterministic real symmetric matrix such that ||F|| < n'^ for 
some 7 > 0. Suppose b is an independent random vector whose coordinates are iid 
copies of let p be a deterministic vector. Then, for any a > 0, there exists C > 0 
(depending only on a, 7, f, and () such that (IT -\- F,b p) is controllable with 
probability at least 1 — Cn~°‘. 

Another interesting problem is the controllability of the pair (IT + F, u) when u 
is uniformly distributed over the unit sphere, independent of IT. 

Theorem 3. Let f be a sub-gaussian random variable with mean zero and unit 
variance; let ( be a (possibly degenerate) sub-gaussian random variable with mean 
zero. Assume W is an n x n Wigner random matrix with atom variables f and (. 
Let F be an nxn deterministic real symmetric matrix such that ||F|| < n'^ for some 
7 > 0. Suppose u is an independent random vector uniformly distributed on the 
unit sphere in K". Then, for any a > 0, there exists C > 0 (depending only on a, 
7, and () such that (IT -\- F,u) is controllable with probability at least 1 — Cn~'^. 

As a consequence of the previous two theorems, we have the following interesting 
result on the controllability of random graphs. 

Corollary 2. Fix 0 < p < 1 and a > 0. Let A be the adjacency matrix of G{n,p). 
Let b be an independent random vector whose entries are iid zero-one Bernoulli 
random variables with parameter p, and let u be an independent random vector 
uniformly distributed on the unit sphere in K.". Then there exists C > 0 (depending 
only on p and a) such that {A, b) and {A, u) are controllable with probability at least 
1 - Gn-°‘. 


3.3. Results for orthogonally invariant ensembles. Due to the continious na¬ 
ture of the orthogonally invariant ensembles and the singular nature of controlla¬ 
bility, we prove stronger universality results for random matrices drawn from this 
ensemble. 

Theorem 4 (Controllability of any orthogonally invariant ensemble). Let M be an 
nxn random real .symmetric matrix drawn from an orthogonally invariant ensemble. 
Let b be any nonzero vector in R.”. Then {M, b) is controllable with probability one. 

As a consequence of the above result, we immediately obtain the following corol¬ 
lary. 
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Corollary 3. Let M he an n x n random real symmetric matrix drawn from an 
orthogonally invariant ensemble. Let b be an independent random vector in R" 
(with any distribution). Then (M, b) is controllable with probability 1 — P (6 = 0). 


4. Tools 


We now present the tools we will need to prove our main results. 


4.1. Tools from linear algebra. Recall that the spectrum of an n x n Hermitian 
matrix A is said to be simple if all of its n eigenvalues are distinct; that is, 

Xi{A) < ■ ■ ■ < XniA). 

In the sequal, we make use of the following criteria for controllability, which is 
based on the Popov-Belevitch-Hautus Controllability Test (see, for instance, [H 
Theorem 12.3]). 

Lemma 1 (Controllability criteria). Let A be an nxn Hermitian matrix, and let b 
be any nonzero vector. If A has simple spectrum and Vj{A)*b ^ 0 for all 1 < j < n, 
then { A, b) is controllable. 

Remark 2. In general, the eigenvectors ui(A),... ,Vn{A) of A are not uniquely 
determined. However, if the spectrum of A is simple, then the eigenvectors of A 
are uniquely determined up to phase. In particular, the condition Vj{A)*b 7 ^ 0 is 
not affected by this choice of phase. 


Proof of Lemma [7J In order to reach a contradiction, assume the matrix 


B-.= {b Ab ■■■ A^-^b) 

does not have full rank. Then there exists a nonzero vector a = (aj)))!^ such that 
Ba = 0. 

For notational convenience, let Ai < • • • < A„ be the eigenvalues of A with 
corresponding orthonormal eigenvectors vi,... ,Vn. By the spectral theorem, for 
any fc S N, 

n 

AH = ^(A^;5)u,. 

i=i 

Thus, using the notation that A^ is the identity matrix, we have 

n—1 n—1 n n /n—1 \ n 

0 = Ba=Y^ OkA^b = ^ ofc ^ H H H > 

fc=o fc=i j=i j=i Vfc=o / i=i 

where 

n—1 

Pr-=T.^]bX]ak. 

k=0 

Since ui,..., are linearly independent, it must be the case that (Ij = 0 for all 
1 < j < We observe that 


n—1 





A, 

VArv 
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Since v*b ^ 0 for all 1 < j < n by assumption, it must be the case that 

( 

a = 0 

V\”“V 

for all 1 < j < n. However, this implies that the matrix 

/A? AO ••• A0\^ 

Ai A 2 • • • A„ 

A . 

Ar^ ArV 

is singular. As A is a Vandermonde matrix, we find that 

det A = ]^(Aj — Ai) = 0. 

i<j 

This implies that the spectrum of A is not simple, a contradiction. □ 


Cauchy’s interlacing theorem implies that the eigenvalues of the Ah minor of 
a Hermitian matrix interlace with the eigenvalues of the original matrix; see [3J 
Corollary III.1.5] for details. One of our main observations is that, if this interlacing 
is strict, then the pair (A, e^) is controllable. For this, we will need the following 
result from |^; see also [HI Lemma 41]. 


Lemma 2 ([Bj; Lemma 41 from HZ])- Let 

An = 


An — 1 
X* 


X 

d 


be an n X n Hermitian matrix, where A„_i is the upper-left (n — 1) x (n— 1) minor 

of An, X € and d £ M.. Let be a unit eigenvector of Xi{An), where 

V € C"“^ and x £ C. Suppose that none of the eigenvalues of A„_i are equal to 
Xi{An). Then 


1 + (Aj(A„_i) - K{An))-^\vj{An-iyx\^ 

where Vj{An-i) is a unit eigenvector corresponding to the eigenvalue Aj(A„_i). In 
particular, under the assumptions above, x ^ 0. 


The above lemma implies that if Xj{An-i) — Ai(A„) ^ 0 for all i,j, then e„ is 
not orthogonal to any eigenvector of A„. The choice of the nth minor (and the nth 
entry of an eigenvector) is arbitrary; indeed, the result can be applied to any minor 
by conjugating the matrix An by a permutation matrix. 

We will use the following lemma in order to show that Aj(A„_i) — Ai(A„) ^ 0 
for all i,j. 




Lemma 3. Let 
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be an nxn Hermitian matrix, where An-i is the upper-left (n— 1) x (n— 1) minor of 
An, X G and d G R. If an eigenvalue of An corresponds with an eigenvalue 

of An-i, then there exists a unit eigenvector w of An-i such that X*w = 0. 


Proof. Assume A is an eigenvalue of A„ with unit eigenvector J , where v G C" ^ 

and q G <C. In addition, assume A is an eigenvalue of An-i with unit eigenvector 
u G From the eigenvalue equation 



we obtain 


An-lV qX = Xv, (3) 

X*v qd = Xq. (4) 

Since = Xu*, we multiply (|31) on the left by u* to obtain qu*X = 0. In 

other words, either u*X = 0 or g = 0. 

If u*X = 0, then the proof is complete (since it is a unit eigenvector of An-i). 
Assume q = 0. Then d is a unit vector, and equations m and dH) imply that 

An-lV = Xv, 

X*v = 0. 

In other words, d is a unit eigenvector of A„_i and X*v = 0. □ 

4.2. Tools from probability theory. We present some probabilistic results we 

will need to prove our main results. The first lemma below is a standard result; we 

include a proof for completeness. 

Lemma 4. Let v G M" be a unit vector, and assume u is a random vector uniformly 
distributed on the unit sphere in M" . Then 

F (v'^u = O) = 0. 

Proof. Observe that for any nxn orthogonal matrix U, v'^u = {Uvf^ {Uu) and Uu 
has the same distribution as u. Thus, it suffices to show that 

P {eju = O) = 0. 

We note that u has the same distribution as 

/ ^n ^ (^1’ ■ • ■ > Cn) ) 

V 2.^fc=i sfc 

where ^i,-.-,^n are iid standard normal random variables. Indeed, this follows 
since the normal vector (^i,..., is rotationally invariant. As J2k=i ^ with 
probability one, we find that 

P(e^"u = 0) = P(Ci = 0) = 0, 

and the proof is complete. □ 

We will also need the following bound for the spectral norm of a Wigner matrix. 
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Lemma 5 (Spectral norm of a Wigner matrix). Let ^ he a sub-gaussian random 
variable with mean zero and unit variance; let be a (possibly degenerate) sub- 
gaussian random variable with mean zero. Assume W is an n x n Wigner random 
matrix with atom variables and (. Then there exists constants Co,co > 0 (de¬ 
pending only on the sub-gaussian moments of and () such that ||W^|| < Cgy/n 
with probability at least 1 — Co exp(—con). 

Proof. It follows from [51 Proposition 5.10] that, for any unit vectors u,v G M" and 
any t > 0, 

¥{\u'^Wv\ >t) < C'exp(-ct2), 

where C,c > 0 depend only on the sub-gaussian moments of ( and (. Thus, by a 
standard net argument (see [T51 Fact 11] for details), we obtain 

P(||1F|| > < C'exp(-c'n) 

for some constants C ,c' > 0 depending only on the sub-gaussian moments of f and 

c. □ 


5. Proofs of the Main Results 

5.1. Proof of results in Section I3.1L This section is devoted to the proof of 
Theorem [T] and Corollary |TJ In view of Lemma [l] it suffices to show that, with 
high probability, the spectrum oi M := W -\- F is simple and Vj{M)'^ei 0 for all 
1 ^ i ^ It was recently shown in m that M has simple spectrum with high 
probability. 

Theorem 5 (Theorem 2.6 from [H]). Let he a sub-gaussian random variable with 
mean zero and unit variance; let ( be a (possibly degenerate) sub-gaussian random 
variable with mean zero. Assume W is an n x n Wigner random matrix with atom 
variables ^ and Q. Let F be an n x n deterministic real symmetric matrix such that 
||P"|| < nA for some 7 > 0. Then, for any a > 0, there exists C > 0 (depending 
only on a, 7 , and () such that W -\- F has simple spectrum with probability at 
least 1 — Cn~°‘. 

Remark 3. [121 Theorem 2.6] is actually more general than what is stated above. 
Indeed, the results in [12] also control the gaps between the eigenvalues X.+i{W + 
F) — \i{W F) for all 1 < j < n — 1. 

We now verify the following result for the eigenvectors oi M ■.= W -\- F. 

Lemma 6 . Let ( be a sub-gaussian random variable with mean zero and unit vari¬ 
ance; let C, be a (possibly degenerate) sub-gaussian random variable with mean zero. 
Assume W is an n x n Wigner random matrix with atom variables ^ and (. Let 
F be an n X n deterministic real symmetric matrix such that jjFjj < n'*' for some 
7 > 0. Then, for any a > 0, there exists C > 0 (depending only on a, 7 , f, and () 
such that, 

P (3z, j S {1,..., n} such that Vj{M)'^ei = O) < Cn~°‘, 
where M :=W -\- F. 

In order to prove LemmalHl we will need the following notation and results. Let 
X = {xk)k=i G R" and let ^ be a real-valued random variable. Consider the small 
ball probability 

Pi,six) := SUpP (l^lXl -I- f,nXn - 

aGR 


a] < (5), 


(5) 
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where are iid copies of In particular, even if the vector x is random, 

we assume ,..., are independent of x and the probability above is taken with 
respect to the random variables In other words, we condition on (or 

freeze) the coordinates of x when we compute the small ball probability in ([5]). 

We will need the following result concerning the small ball probability when x 
happens to be an eigenvector of a Wigner matrix. 

Theorem 6 (Theorem 4.3 from [12]1. Let ^ be a sub-gaussian random variable 
with mean zero and unit variance; let ( be a (possibly degenerate) sub-gaussian 
random variable with mean zero. Assume W is an n x n Wigner random matrix 
with atom variables ^ and ((. Let F be an n x n deterministic real symmetric 
matrix such that ||F|| < n'^ for some 7 > 0. Then, for any a > 0, there exists 
/3,/3o, C”, Co > 0 (depending only on a, 7 , and (() such that the following holds 
with probability at least 1 — Co exp(—/3on).' every unit eigenvector v ofW-\-F obeys 
the anti-concentration estimate 

P^,n-piv) < Cn~°'. 

The proof of Theorem [S] presented in m is based on the inverse Littlewood- 
Offord theory introduced in [18]; see m for a survey. In broad terms, Nguyen, 
Tao, and Vu establish Theorem |6| using the following procedure. It is first shown 
that if p^^s{x) > Cn~°‘, then the vector x must have a rich additive structure. 
Intuitively, this follows since the sum 

n 

can only be concentrated on a given value when most of the coefficients Xk are 
arithmetically well comparable. On the other hand, if we take x to be an eigenvector 
of a Wigner random matrix, then one expects the vector x to look random and not 
have any rigid structure. This explains the intuition behind Theorem |6l 

We also note that Theorem [6] is a slight extension of Theorem 4 . 3 ]. We have 
not included the assumption that the eigenvalues of IT + F are no larger than some 
polynomial in n since this follows immediately from the bound ||IT+F|| < ||IT||+n'*' 
and Lemma [5] 

We now prove Lemma [6] 

Proof of Lemma\^ Set M := W -\- F, and let a > 0. Fix 1 < i, j < n. We will 
show that 

F{vj{Mfei = 0) < Cn““"^. 

The claim will then follow from the union bound over 1 < i, j < n. 

Let M', IT', F' denote the (n — 1) x (n — 1) sub-matrices of M, IT, F formed 
by removing the ith row and ith column. Let X, Y, Z denote the Ith columns of 
M, IT, F with the ith entries removed. In particular, M' = IT'-|-F' and X = Y-\-Z. 
In addition, the entries of Y are iid copies of 

Suppose Vj{Mf^ei = 0. Then, by Lemma [2l it must be the case that one of the 
eigenvalues of M' corresponds to Xj{M). Thus, by Lemma [3] there exists a unit 
eigenvector v of M' such that X'^v = 0. Therefore, it suffices to show that 

P (there exists a unit eigenvector v of M' such that X"^v = O) < . 

Observe that M' and X are independent. In other words, conditioning on M' 
does not affect X. Also, recall the definition of the small ball probability given in 
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Let P, C be positive constants to be chosen later. For a fixed realization of 
M', we say that M' is in a (/3, C')-good configuration if 

for all unit eigenvectors v of M'. Let be the event where the spectrum of M' 
is simple and M' is in a (/3, C')-good configuration. Theorems O and [5] imply that 

there exists P,C > 0 such that P . 

For notational convenience, let S denote the set of unit eigenvectors M'. When 
the eigenvalues of M' are simple, the (unit) eigenvectors of M' are uniquely deter¬ 
mined up to sign. Thus, on the event ^ is a set consisting of n — 1 orthonormal 

vectors (where the choice of sign for each vector is arbitrary). Moreover, for a fixed 
realization of M' and for a fixed vector v £ S, X'^v = 0 implies that \Y*v—a\ < n~^ 
for some a G K. In other words, 

= 0) < P5_„-^(u). 

Thus, by conditioning on we apply the union bound over all vectors in S to 

obtain 

P (3v £ S such that X'^v = O) < P (3v £ S such that X'^v = 0|O/3^c) + P(^^,c) 

< Cnn-’^-^ + Cn-'^-^ 

and the proof of the lemma is complete. □ 


We are now ready to prove Theorem [T] 


Proof of Theorem[l\ Set XI := W + F, and let a > 0. Define the event 

D := {the spectrum of M is simple} ^^{vj{M)"^ei ^ 0 for all 1 < i, j < n}. 

By Theorem [5] and Lemma [ 6 l there exists C > 0 such that P(r2'^) < Cn~°‘. On the 
event fl, it follows from Lemma[T]that (M, ep is controllable for each 1 < i < n. □ 

Corollary [T] now follows from Theorem [TJ 


Proof of Corollaryl^ Let 0 < p < 1, and let A be the adjacency matrix of G{n,p). 
Set tr^ := p(l — p). Note that {A,ei) is controllable if and only if {a~^A,ei) is 
controllabe. 

Define the random variables C •= 0 and 

- _J , with probability p, 

\ — with probability 1 — p. 

In particular, ^ and C both have mean zero and ^ has unit variance. In addition, 
both are sub-gaussian random variables (with sub-gaussian moments depending 
only on p). Let IF be an u x n Wigner random matrix with atom variables ^ and 
C Define the deterministic matrix F = {fij)'f j^i by 


fij 


0 , if i = j. 


Then a~^A has the same distribution as IF + F. Thus, it suffices to estimate the 
probability that (IF + F,ei) is controllable for all 1 < f < n. Hence, the proof is 
now complete by applying Theorem [TJ □ 
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5.2. Proof of results in Section 13.21 This section is devoted to the proof of 
Theorem [21 Theorem |3l and Corollary |2l Unsurprisingly, the proofs presented here 
are very similar to the proofs of the previous section. 

Proof of Theorem[^ Set M := W + F, and let a > 0. In view of Lemma [Hand 
Theorem O it suffices to show that 

P (there exists a unit eigenvector w of M such that + /r) = O) < Cn~°‘. 

We proceed as in the proof of Lemma [H For a fixed realization of M, we say 
that M is in a (/?, C')-good configuration if 

Pi,n-p{y) < 

for all unit eigenvectors v oi M-, here denotes the small ball probability 

defined in ([5]). Let n/3,c be the event where the spectrum of M is simple and M is 
in a (/3, C')-good configuration. Theorems |S] and [S] imply that there exists /3, C > 0 

such that P 

For notational convenience, let S denote the set of unit eigenvectors M. When 
the eigenvalues of M are simple, the eigenvectors of M are uniquely determined 
up to sign. Thus, on the event S' is a set consisting of n — 1 orthonormal 

vectors (where the choice of sign for each vector is arbitrary). Moreover, for a 
fixed realization of M and for a fixed vector v G S^ {b + ^) = 0 implies that 
\v^b — a| < n~^ for some a G K.. In other words, 

P {v^{b + /i) = 0) < p5_„-/3(v). 

Here we have used the assumption that b is independent of M. Thus, by condition¬ 
ing on we apply the union bound over all vectors in S to obtain 

P (du G S with {b + fP) = O) < P (3n G S with (b -I- /i) = -|- P(n^ 

and the proof is complete. □ 

Proof of Theorem\^ Set M ■.= W + F, and let a > 0. By Lemma (U it suffices to 
show that, with probability at least 1 — Cn~°^, the spectrum of M is simple and 
that Vj{M)'^u 7 ^ 0 for all 1 < j < n. The claims follow from Theorem|5]and Lemma 

El □ 

Using Theorems [2] and |3l we now prove Corollary |21 

Proof of Corollary\^ We proceed as in the proof of Corollary [T] Let a > 0. Let 
0 < p < 1, and let A be the adjacency matrix of G{n,p). Assume 6 is a random 
vector whose entries are iid zero-one Bernoulli random variables with parameter p. 
Set (7^ := p(l —p). Observe that (A, 6) is controllable if and only if {a~^A,h) is 
controllable. Thus, it suffices to compute the probability that {a~^A,b) is control¬ 
lable. 

Define the random variables C := 0 and 

^_ J , with probability p, 

^ ' I probability 1 — p. 

In particular, ^ and ^ both have mean zero and ^ has unit variance. In addition, 
both are sub-gaussian random variables (with sub-gaussian moments depending 
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only on p). Let W he an nx n Wigner random matrix, independent of b, with atom 
variables ^ and (. Define the deterministic matrix F = by 


fij 


0 , if i = j. 


Then a~^A has the same distribution as W + F. Thus, it suffices to estimate the 
probability that {W + F, b) is controllable. Since rank is invariant under scaling, it 
in fact suffices to compute the probability that {W + F,a~^b) is controllable. 

Let b' be a random vector, independent of VL, whose entries are iid copies of 
let p he a vector in K" whose entries are all taken to be p/<T. Then, a~^b has the 
same distribution as b' + p. Therefore it suffices to estimate the probability that 
{W + F,b' + p) is controllable. By Theorem [21 we conclude that {W + F,b' + p) is 
controllable with probability at least 1 — Cn~°‘. 

We now show that {A, u) is controllable, where u is a random vector uniformly 
distributed on the unit sphere, independent of A. By the arguments above, it 
suffices to estimate the probability that [W + F,u) is controllable, where W and 
F are defined above and W, u are taken to be independent. The claim now follows 
from Theorem |3l □ 


5.3. Proof of results in Section 13.31 This section is devoted to the proof of 
Theorem [Hand Corollary [31 

Proof of Theorem Let M be an n x n real symmetric matrix drawn from an 
orthogonally invariant ensemble, and let 6 be a nonzero vector in R". Since the 
rank of a matrix is invariant under scaling, it suffices to assume that 6 is a unit 
vector. 

Let 17 be an n X n random orthogonal matrix distributed according to Haar 
measure, independent of M. Since M is orthogonally invariant (see, for example, 
[H Section 2.2]), it follows that UMU* has the same distribution as M. Thus, 

P (rank {b Mb ■■■ M^-^b) = n) 

= P (rank (1717*6 UMU*b ■■■ U M^-^U*b) = n) . 

As U is an invertible matrix, it suffices to show that (M, U*b) is controllable with 
probability one. 

Since 17*6 is a random vector uniformly distributed on the unit sphere in R", it 
in fact suffices to show that (M, u) is controllable with probability one, where u is 
uniformly distributed on the unit sphere, independent of M. 

It is well-known (see, for instance, [U Section 2.3]), that the spectrum of M 
is simple with probability one. By Lemma (T] it thus suffices to show that, with 
probability one, Vj{M)'^u 0 for all 1 < j < n. The claim now follows from 
Lemma [U □ 


Proof of Corollary\^ Let Tl be the event that (M, 6) is controllable. We condition 
on 6 (since M and 6 are independent) to obtain 

P (D) = P (D|6 ^ 0) P(6 0) -b P (D|6 = 0) P(6 = 0) 

= P(6 7^ 0) 

= 1-P(6 = 0) 
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by Theorem U) Here we used the fact that {A, 0) is not controllable for any matrix 
A. □ 
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6. Closing remarks 

In many cases, the above results can also be extended to Hermitian random ma¬ 
trices whose off-diagonal entries are complex-valued random variables. Although, 
the focus of this work has been the study of SISO systems, some generalizations to 
the case of MIMO systems is fairly straightforward. In addition, although our study 
is focused on controllability of random systems, due to the duality of controllability 
and observability, all of our results can be formulated naturally for observability 
and controllability of random systems. 

To the best of the authors’ knowledge, this work is the first study of control¬ 
lability of random systems and many questions are left for future works. One of 
the main open problems is to understand the extent to which the universality of 
controllability and minimal controllability holds. For example, the independence 
assumption on the entries of the Wigner matrices discussed in this work is essen¬ 
tial to our analysis. This is certainly not the case when one considers Laplacian 
dynamics. The existence of such universality results for other ensembles of random 
matrices is left open for future studies. Also, controllability of a random matrix 
with an arbitrary deterministic input vector other than (such as the case of 
Conjecture [2]) remains unanswered. 
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